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Abstract 

First, we show that implied normal volatility is intimately linked with the incomplete 
Gamma function. Then, we deduce an expansion on implied normal volatility in terms 
of the time-value of a European call option. Then, we formulate an equivalence between 
the implied normal volatility and the lognormal implied volatility with any strike and any 
model. This generalizes a known result for the SABR model. Finally, we adress the issue 
of the "breakeven move" of a delta-hedged portfolio. 



JEL Classification G12 G13 C65 
Keywords : Smile asymptotics, implied normal volatility, "breakeven move" . 

1 Introduction 



This article comes within the scope of the study of the asymptotics of implied volatility 
which has been considered extensively ([E], [11], [?])• Asymptotics of implied volatility 
are important for different reasons. First, they give information on the behaviour of the 
underlying through the moment formula ([H]) or the tail-wing formula ([4J). Second, they 
allow a full correspondence between vanilla prices and implied volatilities. With such a 
correspondence, asymptotics in call prices can be easily transformed into asymptotics in 
implied volatilities. When applied to a specific model, asymptotics are widely used as smile 
generators (|13j). In practice, other models are then used for pricing options using tools like 
Monte- Carlo simulations. 

So far, all the asymptotics studied by authors concern asymptotics for implied lognormal 
volatility. In this article, we consider implied normal volatility which refers to the Bachelier 
model. Why is it interesting to consider normal implied volatility? First, for short maturities, 
the Bachelier process makes more sense than the Black-Scholes model. Indeed, the behaviour 
of the underlying from one day to another is generally well approximated by a Gaussian 
random variable (see [2Q])- That's the reason why the Bachelier model is very popular in high 
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frequency trading ([3]). Second, the "breakeven move" of a delta- hedged portfolio option 
is easily interpreted as normal volatility. Generally, the P & L of a book of delta-hedged 
options is positive if the (historical) volatility of the underlying is greater than a breakeven 
volatility which has to be expressed in normal volatility. Moreover, it makes more sense to 
compare implied normal volatilities with historical moves of the underlying as can be done 
by a market risk department. Likewise, some markets such as fixed-income markets with 
products like spread-options are quoted in terms of implied normal volatility (|15j). Finally, 
the skewness of swaption prices is much reduced if priced in terms of normal volatility instead 
of lognormal volatility. Therefore, it is important to have a robust and quick way to compute 
implied normal volatilities from market prices and also to be able to switch between lognormal 
volatilities and normal volatilities. 

What kind of asymptotics should we consider? Most of the approximations in option 
pricing theory are made under the assumption that the maturity is either small (see the 
Hagan et al formula |13j for instance) or large (|10j); it is actually assumed that a certain 
time-variance a 2 T is either small or large. A possible way to derive such approximations 
is to replace the factor of volatility a by ea and then set e = 1. This can be done at the 
partial differential equation level (see all the techniques coming from physics |13| ) as well as 
directly at the stochastic differential equation level with the help of the Wiener chaos theory 
for instance |21j. Other types of asymptotics are obtained by considering large strikes. In 
our approach, we unify all these types of asymptotics (see JTT] and [12] for the lognormal 
case). Indeed, we obtain an approximation of the implied normal volatility as an asymptotic 
expansion in a parameter A for A<1 and it turns out that A — > when T — > or K — > +oo. 

This study is organized as follows. We first give another expression for the pricing of a 
European call option which involves an incomplete Gamma function (Proposition [T]). Then, 
we inverse this function asymptotically and obtain an expansion of normal implied volatility. 
This is particularly important if we want to quickly obtain the implied normal volatilities 
from call prices as is the case in high frequency trading ([3]). The formula is also potentially 
useful theoretically if, given an approximation for the price of a European call option or a 
spread option (for instance in the framework of the Heston or the S ABR model) , we want to 
obtain an approximation of the normal implied volatility. Finally, we restrict our formula to 
the order and we compare it to a similar formula for the lognormal case. Then, we obtain 
an equivalence between normal volatility and lognormal volatility. We use it also to compare 
the Black-Scholes greeks to the Bachelier greeks. Finally, we consider a delta-hedged portfolio 
and we compute the "breakeven move" in the normal case as well as in the lognormal case. 

2 Another pricing formula for call options in the Bachelier 
model 

In a Bachelier model, the dynamic of a stock (St) is given by: 

dS t = a N dW t , 

with initial value S at t = 0. The so-called normal volatility is related to the price of 
a call C (T, K ) struck at K with maturity T by the formula [20J : 



C(T,K) = (S-K)N 
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Following Ropper-Rutkowski f[19j). we can isolate the volatility a at in the pricing formula. 

Definition 1 Let us denote by TY(K,T) (or simply TY) the time-value of a European call 
option struck at strike K with maturity T: TV(T, K) := C(T,K) — (S — K) + . 

Proposition 1 In the Bachelier model, 

(-i,&Sn if k + s 
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TV(T,K) = { ^ V 2 ' 2^ry ^ ' (2) 
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where T (a, z) is i/ie incomplete Gamma function: 



otherwise 



r+oo 

T(a,z)= / n a_1 exp(— n)dn 
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The proof is given in the Appendix. It is clear from Proposition [T] that the real- valued 
function T i-> C(T,K) is non-decreasing, positive, C(0, -K") = (£? — iO+ and lim C(T,K) = 

T— s>+oo 

+oo. So, given the price of a European call option C, there is a unique real number ctn(T, K) 
such that C(T,K) = C with a normal volatility <tjv = cf^{T,K). We say that o~n(T,K) is 
the normal implied volatility. 

Note 1 We remark that — — — depends only of — — — (only one variable). 

One of the interests of Proposition Q] is that there are efficient algorithms to compute 
the inverse of the incomplete Gamma function. In particular, it is implemented in Matlab. 
Therefore, it is always easy to get the implied normal volatility from call prices ([16J) . Such 
a task is not always easy in the lognormal case (|14j). especially when we are far from the 
money. 

Corollary 1 Let p be an integer. Then, 



with 



(2p + l)\f a 2 N T 



p\ \(S-K) 2 

This last equation comes naturally from a well known asymptotic expansion of T(a, z) for 
large z (see Formula 6.5.32 in ([I]))- 
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Note 2 From either pricing formula (OP or we can notice that we can use the same 
demarch to price large strikes and short maturities European options (as expansions in both 
cases are similar)... This comes from the fact that: 

C(X 2 T, XS + (1 - A) if) = AC(T, K) 
for any non-negative real A. This is particular to the Bachelier model. 

To compare with the lognormal case, it can be advantageous to introduce the following 
notations. 

(Tat VT 

Definition 2 For K ^ S, we set (the symbol N stands for "Normal") On := — — — , xn '■= 

K . . fA^\ 29% (A^\ . 1 
. um := — 7m := In 



— - - 1, 7jv := In - ) ,u N := -f N := In , A :— 

S \\xn\J x% \\xn\J \n( Ty( l' K) 

Then, by Corollary [Q for K ^ S and p £ W, 

^i^ = 4.-*g^ a i^ +J #) (3) 

with R f '} e O (<ft\ and a?' - (2< = + 1)1 
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In comparison to the normal case, in the lognormal case, we have (|12|): 



40Fe-^TV(r,if) _ 3 ^ (_!)* (fe) (p) \ 

with u LA r := &LN ■= (JlnVt, x ln := In y— J , i?( iAr) G O (#1^) , 

and (2Jfe + 1)!! := JJ(2j + 1). Here, ctln denotes the lognormal implied volatility and "LN" 
stands for "Lognormal". 

3 Asymptotics of the implied normal volatility 

Let us assume that K ^ S. Using ©, we get: 



e 4 (j2 a k u k N + O {u p N ) \ = e™ e~* 
\fc=o / 
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with ak '■= — r — afc. Therefore, by Lemma 1 of 1121 . we get the following proposition. 
2 
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Proposition 2 Let us denote by TV the time-value of a European call option, o~n its implied 
normal volatility and T the maturity of the option. Set A := tv^' ^ n := ^ n ' 



and xn 



K 



\\ x n\J 



1. Let us assume that K ^ S. Then, in the case when T — > 0, we have the 



(S — K) 2 

following expansion for the time-variance of the call option: o~ N T = un with 



u N = \-h 2 ln\ + 1N X 2 + ~ 4 \ 3 ln 2 (A) + (j - 3 7 tv) A 3 ln(A) + ^ - *- lN + A 3 + o (A 3 ) 

(6) 

In the lognormal case [12], for short expiries, the asymptotic expansion of o~ 2 LN T is given by 

ln 2 (4) 
crf^T = u LN with 



uln = A-^A 2 lnA+ 7L7V A 2 +^A 3 ln 2 (A)+^ - 3 7L ^ A 3 ln(A)+ (j 2 LN - ^ lLN - a'^j A 3 +o (A 3 ) 
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First, we note that A = uln + o(uln)- Then, comparing the two results © and ([7]) for 
K ^ S, we obtain: 



UN = U L N + {iN ~ ILn) u\ n + O (u LN ln(li LA r)) 



So, 
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Since = — > 0, we deduce that 

xln hi 6 — hi K 



(8) 



Moreover, we have: 



Hence, we get: 
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Note also that at the money, the situation is quite easy. On the one hand, we have (cf 
Proposition [I]) 



ON 




(11) 



On the other hand, we have ([12J, Proposition 2): 

2V2 , 



C = erf 



Therefore, we state the following result. 
Corollary 2 • If K 7^ S , we have: 
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In particular, <jn ~ Sctln when T — > 0. 



Note 3 When K — >• S", u>e can check that 




In other terms, for K ^ S, 
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where m 



A' 



xtv + 1 is the moneyness. 



This formula was known (even if it was not stated explicitly) in the SABR model (see 
the Hagan et al formula |13j). By differentiating Formula ([3]) with respect to m, it turns out 
that the Black-Scholes skew ^Sf" at the money (m = 1) generated by the Bachelier model 



1S d< drn = ~11§ L ( a LN is by definition the implied lognormal volatility). Therefore, the 
Bachelier model is highly skewed ATM (a slope of —50% x ^-)- Another way to explain this 
feature is that given call prices, when we use the BS model, the function <jln is a decreasing 
and convex function of m, i.e., it generates a skew, while the function is a rather flat 
function of m. Thus, normal volatility is most suited for products such as swaptions for 
instance. 



4 Comparing greeks and delta-hedged portfolios 

Let us denote by Ajy, Tjy, vjy, On (resp. Aln, ^ln, ^ln , ®ln), the delta, gamma, vega and 
theta in the Bachelier (resp. Black-Scholes) model. For instance, vn = By differentiating 
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(HD, we get: 

A n = n(^L\ (10, 



On the other hand, it is known that: 



liv = .v|!22z** + i VT] (10) 



So, by Corollary [51 we get: A/jv ~ A^r for a maturity T<1. By differentiating we 
obtain: 

a N VT \a N VTj 
In the Black-Scholes model, we have:rrrrrr 



1 {\nS-lnK 1 /-\ 

Hence, with the help of Corollary [21 
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Now, we consider the Vega. It is shown in the Appendix (Formula (JH)) that: 

C(T, K) = (S- K)+ + S n du (10) 



o 
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So, differentiating by ctn, we get: 



VN = S x l^e W ( ( i 2 K f ) (10) 
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In contrast, the vega in the Black-Scholes model is 



T ( OnS-\nK + ±aj N Ty 

2 - 
LN J 



In the same way, we can compare the two thetas. So, we get the following proposition. 
Proposition 3 When T — > 0, and under the hypothesis of bounded volatilities, we have: 
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The first equivalence A at ~ Aln shows that hedging in the Bachelier framework is more 
or less like hedging in a Black-Scholes framework. However, the "breakeven move" of a delta- 
hedged portfolio is not the same. By definition, the "breakeven move" of a delta-hedged 
portfolio is the number [i such that over a short horizon St, P&L > if the change in S is 
> /U. In general, we have: 



P&L = -est + ^T(AS) 2 (15) 
= It [(AS) 2 - /x 2 ] (16) 



So, with St = 1, 



/29 , x 

» = yjr- (16) 

Using Proposition [31 we get that the "breakeven move" fiLN in the Black-Scholes model 
is related with the "breakeven move" ///v in the Bachelier model by: 

lnm 

[i-LN = tUn (16) 

m — 1 

with m = 7^. So, at the money, /ilat ~ ^at. However, if K < S (resp. K > S) then uln > £W 

ln(?Ti) 

(resp. un > Uln)- We have represented in Fig 1, the function m \-t which gives the 

m — 1 

smile of the Bachelier model (cf Corollary [2]) as well as the ratio of "breakeven moves" LN . 

UN 

So, depending on the view of the trader on the short term dynamic of the underlying 
(normal or lognormal diffusion), he will adjust or not the "breakeven move" of its delta- 
hedged portfolio by the factor -^j- 
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Figure 1: The smile of volatility of the Bachelier model (normalized by the factor where 
cjat denotes the normal volatility) and the ratio of "breakeven moves" uln/hn- 
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Appendix 

Proof of Proposition [TJ. We have: 
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We have: 



^-fo") : = — — n ( — ) +n ( — ) +n ( — ^-)-( — — ) n 



= n 



W \ u / \ u I \ u * I \ u I \ u 
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We have used: n'(£) = -£n(£). Since /(£,0) = = (1 - £)+, w e deduce that: 



Set 



Then, 



/(«,») = (i-0 + + jf»(^) 



C ( r - K )=(l- e)+ + F^ ^ 



Let us assume that both 0/0 and £ / 1. With the change of variable v = we get 



*<«■») = n -«i /J 
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So, with the new change of variable u = ^v 2 , we have: 



i r +oc s 



1 H-Cr^ 1 I 1 "*' 2 



Sl V 2' 20 2 

where T(a, z) is the incomplete Gamma function. At the money, we have simply: 

C(T,K) = aNVT 
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